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Abstract 

We calculate new contributions to the — > tt^uu decay in models where 
neutrino Majorana masses require an extension of the scalar sector. First, we 
study a model where the neutrino mass is induced by the vacuum expectation 
value of an SU (2)-triplet scalar. Second, we study the Zee model where the 
Majorana mass comes from one loop diagrams involving a singly charged, 
5C/(2)-singlet scalar. In both models, the Yukawa couplings that involve the 
new scalar and the neutrinos could be of order one. We find, however, that the 
contributions to the Kl — > tt^uu decay mediated by the new scalars depend 
on the neutrino masses rather than the Yukawa couplings and are, therefore, 
negligibly small. 
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I. INTRODUCTION 



Within the standard model (SM) the Kl tt^ui? decay is known to be a CP violating 
(CPV) process to a very good approximation |l]J^, and subject to a clean theoretical inter- 
pretation 0. In the SM, CP conserving (CPC) contributions are chirally suppressed and 
smaller by many orders of magnitude than the CPV ones 0]. In a previous work we 
calculated the CPC contributions that arise when the SM Lagrangian is extended to include 
neutrino mass terms. These contributions are known to be chirally enhanced, but we found 
that they are suppressed by a factor of order m^/m^ which makes them negligibly small. 
In this work we study models where, in addition to neutrino mass terms, there are new light 
scalars. The question that we ask is whether the fact that such scalars can have Yukawa 
couplings of order one to neutrinos allows for a situation where the new contributions are 
chirally enhanced while avoiding the m^/m^ suppression factor. 

Specifically, we consider the following two models: 

1. Neutrino masses are related to the vacuum expectation value of an S'f/(2)-triplet 
scalar. We present the model and calculate the new contributions to the Kl tx^vv decay 
in section II. 

2. Neutrino masses are related to loop diagrams involving a singly charged SU (2)- 
singlet and an additional S'f/(2)-doublet scalar. We present the model and calculate the 
new contributions to the Kl —>■ tt^uu decay in section III. 

A summary of our conclusions is given in section IV. 

II. AN SU{2) TRIPLET SCALAR 
A. The Model 

We consider the SM Lagrangian with the addition of an SU(2) triplet Higgs field (A). 
The most general scalar potential is given by 

V{H, A) = ^{HlWf + ^(Aj A*)2 + XHjWAlA' + \W^a%WA'H%A" 
2 2 

+ m{A^H'^H'^ + V2A- H+H^ + A—H+H+) + f4jHjH' + /i^Af A* + h.c , (2.1) 

where A* = [A++, A+, A°], W = and a,t are are SU(2) generators in the 

J = 1/2, 1 representations, respectively. The m-term in (|2.1|) breaks the global lepton 
number symmetry (L) so that the phenomenologically unacceptable [^,0 Majoron is 
avoided. For simplicity, we assume that the couplings in ( p.l|) conserve CP. As concerns the 
VEVs, = ^ and (A°) = we assume for simplicity that they are both real and we 
take into account the constraints from the p parameter 

^<10-^ (2.2) 

Vh 
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In order to calculate the contribution to the decay, the scalar mass eigenstates and mixing 
angles should be identified. With the assumption that the scalar potential V{H, A) ( |2.1| ) is 
real, the imaginary and the real parts of the neutral scalars remain unmixed. For the CP 
even fields, \/2TZe{H^) and ^/2TZe{A^), we have the following mass matrix (we neglect 
terms that are higher order in v^/vh and assume that the A's are positive and of order one): 



7^ 



V2 

VH 



A 
'V2 



(A% + m) 



m 



(2.3) 



where A = A + A. The eigenvectors of are: 

^1= ne{H^) COS7 + 7^e(A°) sin7 
^2= -ne{H^) sin 7 + 7^e(A°) C0S7 ^ 

where 



(2.4) 



tan 7 




m 



V2VA \ 



Ah + 



A + 



m 



va 



< 1 



(2.5) 



In deriving the inequality in eq. ( |2.5| ) and below we assume that there are no fine-tuned 
cancellations among the independent parameters of the scalar sector. The eigenvectors 
correspond to the following mass eigenvalues: 



'?1,2 



m ^ 

V2VA \ 



Aff + 



m 



V2VA, 



+ 4- 



A + 



m 



V2rr, 



m 



(2.6) 



>^H + ^' V2VA 

For the CP odd fields, \/2X'm{H'^) and A/22rm(A°), we have the following mass matrix: 

M| = y2m("^J^ . (2.7) 

The eigenvectors of M| are: 



G° = ImiE^) COST] + Jm(A°) sin 77 
J° = — Xm(if°) sin?7 + Xm(A'^) cos 77 . 



(2.8) 



where 
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tanr^l 



2VA 

vh 



< 1. 



(2.9) 



We learn that mixing in the CP-odd sector is very small and we neglect it from here on. 
The eigenvectors correspond to the following mass eigenvalues ||^: 



rriQO JO 



0, 



mvfj 

V2VA 



(2.10) 



where the massless component ~ \/2Xm{H^) corresponds to the unphysical Goldstone 
boson which is eaten by the Z boson, while J° ~ \/2Xm{/\P) corresponds to the would- 
be Majoron. As expected its squared mass is proportional to the explicit lepton number 
violating parameter m. 

For the singly charged scalars, H"^ and A^, we have the following mass matrix: 



Ml 



71^ 71 



-2^ V2 

VH 



(2.11) 



The eigenvectors of are: 



G^= cos 7]' + sin r]' 
= —H^ sin T]' + cos rj' . 



where 



tan rj'l 



V2VA 



Vh 



< 1. 



(2.12) 



(2.13) 



We learn that mixing in the charged sector is very small and we neglect it from here on. 
The eigenvectors correspond to the following mass eigenvalues: 



0, 



X 



(2.14) 



m 

where the massless component ~ corresponds to the unphysical charged Goldstone 
boson which is eaten by the bosons, while ~ A^ corresponds to a new physical 
charged scalar. 

Neutrino Majorana masses are induced via the following interaction terms: 



(2.15) 



where r is an SU(2) generator in the J = 1/2 spinor representation. For simplicity, we take 
fmn = fm^mn with fm real. In this way, we avoid unnecessary complications related to flavor 
mixing and CP violation in the neutrino sector (see discussion in ref. [Q). The interaction 
term ( p.l5| ) induces neutrino masses: 



Maj\ 



(2.16) 



with 2 = 1,2,3. In addition, it generates new contributions to the Kl n uu decay. 
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B. The Kl TT^uy Decay Rate 



The dominant new contributions to the Kl — > t^^vv decay come from the diagrams 
presented in fig. |I|, which generate the following effective Hamiltonian: 



{dsM'rr'v,){xt^+Xr+Xr) 



.(b) 



(2.17) 



where the X^^*'' function corresponds to the diagram in fig 0(1). In the X'f^ and xf^ terms, 
the L breaking effects enter through, respectively, the AVrPF and AHH couplings and is 
related to the spontaneous breaking ((A°) 7^ 0): 



X, 



(a) 



f sr^ . f sin^ 7 cos^ 7 \ / 1 A — A ' 

-m^nis Xi I — ^ \ — I I - + 2xi — — 



X 



(X, - 1)3 



^1 

2xi ln{xi 



m 



€2 



l-x? 



9 



(2.18) 



where Xi = {rrii/MwY and Aj = V*gVid- Since the top quark contribution is dominant, (p.l8|) 
can be simplified: 



m 



^2 



X Xt 



'I + 2xt^\ {2xtHx,) + I - x^) ^ - ^ 



{xt - If 



9^ 



(2.19) 



In the xf^ term, the L breaking effect enters through the A — H mixing and is related 
to both the soft breaking (m 7^ 0) and the nonzero VEV of A*^. The calculation of xj,'^^ is 
simplified by the use of the sdH effective coupling TgdH, [represented by a square in fig. |l](c)] 
which was calculated within the SM in ref. . For our model TgdH should be expressed in 
terms of the appropriate masses and mixing angles. The mixing of the charged scalar can, 
however, be neglected [see eq. ( p.l3|) ] and therefore we use directly the SM calculation: 



sdH 



-A, 



mlms fS 4A 



1287r2 \2 



7: + 



H 



9' 



/2(X,) (I+7') 



where 



f2ix) 



X 



2(1 -x)^ 



X , 2 , 13 

In a; H In a; 

1 — X 1 — X 2 2x 



(2.20) 



(2.21) 



Then X^'^^ is given by 



(c) , I sm2'jVH 
-= -Xtm^ms 



8 wa Vm^^ 



2 ? + ^/2(a;t) 



X2. 



9^ 



(2.22) 
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Note that since the lepton and quark operators in the effective Hamiltonian in 
eq. ( |2.17| ) are scalar operators, the new contributions are CPC 

We are interested in finding the upper bound on the new contributions. We therefore 
focus on the region in parameter space that maximizes them. For all scalar masses, we will 
use a lower bound of 45 GeV, thus avoiding any conflict with constraints from the invisible 
width of the Z boson. From eqs. ( |2.1U| ) and ( ^.141 ) we then find: 

sign(m/t;A) = -1, \m/vA\ ~ 1- (2.23) 



Substituting the proper values for the masses and angles (|2.5|) , (|2.6| ) into the functions 
X^("\xf) (IZISl ) and Xf^ (^) we find: 



x\ ^ + ^ ~ [^^^ v^/'^fi 

Utt \ III 



H 

2 



■'^r'-^^^U-£t^^Oi.l/.l) (2.24) 



mims / Vf 

'^H \^hVh +mvA 

For m ^ vh X^"'^ + X^''^ and X^"^-* are all highly suppressed, together with the total rate. 
Larger contributions are found with m in the intermediate regime ~ ~ "W//, where 



Xi"^ + XP and Xf are all of the order of ^ ~ ^ 



The p^Prn^ coupling, that is fm of eq. (|2.15| ), is proportional to m^^jvA- Naively, one 



may think that this mechanism of inducing neutrino masses can induce a contribution to the 
Ki -nvv rate that is enhanced by a factor of order vh/va compared to the mechanisms 
of 0. We find that this is not the case. For the diagrams in fig. |l|(a,b), there is a va- 
factor in the lyiyA and HHA couplings. For the diagram in fig. |l|(c), there is a factor 
of sin 27t>//(l /m|^ — l/m^J ~ va/v'^h- either case, the final result is proportional to 
fmVA ~ i^um ci'iid there is no enhancement. 

We are now in a position to compare the contribution of the triplet scalar to the leading, 
CPV one 



Ta{,Kl-^ 'n^w) . ( Mk m,, V .^^nf m„ 1^ Ah GeV 



10 MeV 



4 



(2.25) 



ight 

where m^^. stands for the lighter between m^^ and m^^. While the direct bound is rrii,^ < 



18.2 MeV [0], there is a significantly stronger bound from cosmology, rrii, ^ 10 eV [12 
Therefore, very likely, R^py < lO^^s. 



III. THE ZEE MODEL 
A. The Model 

The Zee model enlarges the SM scalar sector by a charged SU{2) singlet 0"'" and an 



SU{2) doublet H2. The SM doublet is denoted by Hi. The Lagrangian of the model is [Il3| : 
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(3.1) 



where L is a lepton doublet, a, (3 are SU(2) indices, m,n are flavor indices, C is the Dirac 
charge conjugation matrix and V{Hi, H2,4>'^) is the most general scalar potential which 
respects the SM gauge symmetries and L [0]: 



.2\ 2 



.2\ 2 



2 



+ A7 



+ A. 



\Hi\' + Xc 



H2\" + M{H^^iT^H'^(f)- + h.c). (3.2) 



Both Higgs doublets develop nonzero VEVs: 



fl,2 



(3.3) 



For = 0, the parameters ai^2 would be equal to f 1,2 respectively. With a nonzero fi term, 
f 1 2 are complicated functions of Oj, Aj and /i. The L-charges of the scalars fields are: 



L{H^) = 0, L{H2) = 2, L{ 



-2. 



(3.4) 



Then we see that L is broken spontaneously via the nonzero VEV of H2. The /i-term in 
eq. ( |3.1|) breaks L explicitly so that the phenomenologically unacceptable Majoron |l7|, p^Jl5 



is avoided. For simplicity, we assume that CP is conserved in the lepton sector and take all 
the dimensionless scalar couplings to be of order one. 

In order to calculate the new contribution to the decay, the scalar mass eigenstates 
and mixing angles should be identified. There are seven physical and three unphysical 
combinations (eaten by the Z and the bosons). Out of the seven physical combinations, 
four are charged and three are neutral. The neutral ones are the would-be Majoron [which 
gains mass due to the fi term in eq. (|3.1|) ] and two other real fields. As long as fi is 
real, the real and imaginary parts of if°2 remain unmixed (as in the cases in which L is 
spontaneously broken |T^|T^), thus the imaginary physical combination corresponds to the 
would-be Majoron which is irrelevant to our calculation. 

The real neutral mass matrix is read from the scalar potential of eq. ( p.ll) : 



7^ 



2Xi{3v{' 



- ai2) + X3{v2^ 
2X3V1V2 + jJ, 



02^) 2A3f if 2 + /i 

2X2{3v2^ - a2^) + Xsivi' 



ai 



(3.5) 



The eigenvectors of M-ji are: 



pi = -JZeiH^) cose + ne{H^) sm9 
P2 = -ne{Hl) sin + 7^e(i^°) cos 6 . 



(3.6) 
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Each eigenvector corresponds to an eigenvalue denoted by rup^ ^. The masses 2 and tan^ 
are comphcated functions of Oj, Aj and /x. 
The charged mass matrix is: 



M, 



c 



\v1V2 + 2^1 \vf - 2fx^ -V2Mvi 
V V2MV2 -V2Mvi 2X6m^ + Xavl + Xgv^ 



(3.7) 



where A = A4 + A5. The massless combination of the charged fields (see e.g refs. [0,|T^) is: 

= Hfcos5 + HtsmS, (3.8) 

with 

(3.9) 



tan S = — . 



The massless combination is not affected by the addition of the singlet field 0^, nor by the n 
term, since it is determined (according to the Goldstone theorem) only by the broken SU{2) l 
generators. The physical charged fields are given by the following linear combinations of 
H^2 0^ (which must be orthogonal to G^): 



-Hf sin 5 + Hf cos 5) cos /? + 0^ sin /3 
[~Hf sin 5 + Hf cos 5) sin /? + 0^ cos (3 . 



(3.10) 



with tan/5 being a complicated function of the scalar potential parameters. The new in- 
teractions in eq. ( |3.1| ) induced a neutrino Majorana mass matrix M^^^. For simplicity we 
assume the dominance of the one loop induced mass |ll3|JI^ , [T6|Jl9(| : 



2y^ 

(4^ 



/m£ tan 5 sin 2j3 



.nip 



mz 



In 



m 



X2 



w 



m 



(3.11) 



XI . 



with m, i flavor indices. 



B. The Kl TT^vu Decay Rate 

The new interactions also generate new contributions to the Kl tt^uu decay. Note that 
the one loop induced mass and the new contributions to the decay (shown in fig. |]) are related 
to mixing between the charged scalars and therefore vanish in the limit 5, /3 — > 0. Below 
we concentrate on the contributions which are dominant when sin 6 is small. Contributions 
that are of higher order do not modify significantly the results, even with sin 5 ~ 1, and 
thus they are omitted. Note that, since vi induces the top mass, we always have tan 5 < 1. 

The scalar operator is induced by the neutral Higgs mediated penguin diagram shown 
in fig. H . The square in the figures represents an effective sdHi^2 vertex denoted by VsdHi^2-i 
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similar to the one we encountered in section II. B . The effective vertices TgdHi 2 induced 
by the Zee model fields differ from the one calculated within the SM |]T^. Since, however, 
we are only interested in finding an upper bound on the decay rate, we can simply set the 
charged mixing angles factors to one, and replace the boson masses in the propagator (as 
explained in section II. B) with m^^^.^^^ ~ 45 GeV. Then we get: 



sdH, 



~ 1287r2m3 

Xlight 



(3.12) 



where CsdH, is a constant of 0(10). Neglecting subdominant contributions of 0[{me/mp^^^.y], 
the diagrams in fig. |^ generate the following CPC effective Hamiltonian: 



n 



y/2 2ti sin^ Bvk 



£m 



^ em J ^ 



(3.13) 



where the Y^^ function corresponds to the diagram in fig. |^(i). Y^^' is given by: 



(a) 



Y, 



(a) _ r^M^mMw m^t^^s 



COS 5 ml 

Xlight L 



C. 



cos^ 6 sin 6 



sdHi 



C. 



sin 2^ 



m 



pi 



m 



sdH2' 



P2 



m 



pi 



m 



P2 . 



(3.14) 



The calculation of the diagram in fig. ^(b) is much more involved. This is due to the fact 
that each of the trilinear couplings HfH^H^ and H^H^cf)'^ k = 1,2) translates into a 
set of eight couplings in the mass basis. In order to estimate the upper bound on the rate we 
set (again) all mixing factors to unity and replace the charged scalar propagators by their 



maximal values, i.e. -p^; 
then find: 



l2 



fc2- 



^2 



fc2- 



l2 

flight 



and 



fc2- 



+ 



^X2 



A;2- 



^2 

flight 



We 



^647r2' 

1 cot 5 
% sin 2/9 



pi Xlight 



M)] 



In^ 



w 



m 



Xi . 



m 



mi 



pi Xlight 

X [Fi {9, mp^ , m^^ , M) + F2 (e, m^, , m^, , M)] 



(3.15) 



where 



Fi(^,mp,,m^,,M) 



/ Til 

CsdH, cos' e + sin' e- 



m 



P2 ■ 



C. 



sin 26 



sdH2' 



2A1 + -A + 2A3 



2A8 + 3 



M 



-1 + 



m 



pi 



m 



P2 . 



(3.16) 



and 



F2{9,mp^,m^^,M) 



I Tfl 

CsdH2 sin'^ + cos'^- 



-C. 



sin 2^ 



m 



sdHi- 



P2 . 



2A2 + 2A3 



1- M 
A + 2A9-5— 

2 Vl 



-1 + 



m 



pi 



m 



P2 , 



(3.17) 
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The M-term in the scalar potential ( p.2| ) does not break any symmetry. One might think 
then that the ratio ^ that appears in eqs. can be arbitrarily large and enhance the 

decay rate. This is however not the case. For very large M, the scalar potential would 
spontaneously break U{1)eai- An even stronger upper bound on M comes from two loop 
contributions to the neutrino masses. We assumed, for simplicity, that the neutrino masses 
are dominated by the one loop contributions. This assumption requires that ^ is not much 
larger than the various A couplings. Adding ( p.l4|) to ( p.l5|) , and applying the approximation 
described above eq. (|3.15|) we find: 



yio-) _|_ yW < 



nil,. M\y mfm 



m 



plight 



'Klight 



C. 



sdH, 



4 + 3- 



w 



C. 



m 



sdH2 



1 + 3- 



Klight 



m 



Klight 



< ^q2 ^vsMw minis 



m 



plight 



(3.18) 



Klight 



where is the largest eigenvalue of M'^^. 

We can now compare the CPC rate of the Zee model with the leading CPV one |TT 

- / ^ ^ 



Zee 
CPV 



Zee I 



CPV 



< 



m 



mz 

(-'light Klight 



10 



-5 



'45 GeV^ 



'45 GeV' 



10 MeV 



m 



plight 



m 



Klight 



(3.19) 



Since, very likely, nii,^ < 10 eV ||20[ , we expect Rcpv ~ 10^^''. Recall that this upper bound 
was obtained using some crude approximations and is expected to be even smaller for an 
exact calculation. 



IV. FINAL CONCLUSIONS 

In this work we examined the question of whether the SM CP violating contribution to 
the Kl — > n^uu decay is still dominant in the presence of new scalars that induce Majorana 
masses for neutrinos. We found the following unambiguous answers: 

(i) For CPC contributions induced by an S'f/(2)-triplet scalar, we get: 



2 



SM 



fr, Pf^ ^10"''- (4-1) 



(ii) For CPC contributions that are generated in the Zee model, we get: 

rZee(^L TtVz/) ^ ( _^ ^^ TTl^Mw mjmK \ < . n-5 
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In obtaining the final bounds in eqs. ( [4.1| ) and {^^ ) we used the direct upper bound on 
m^^ of order 10 MeV. If we use the cosmological bound, nii, < 10 eV, then the bounds 
become stronger by twelve orders of magnitude. It is clear then that the Kl —>■ tt^uu decay 
provides a very clean measurement of fundamental, CP violating properties and that it does 
not probe neutrino masses. 
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